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ABSTRACT. We present a class of graviton-dilaton models which 
leads to a singularity free evolution of the universe. We study 
the evolution of a homogeneous isotropic universe. We follow an 
approach which enables us to analyse the evolution and obtain its 
generic features even in the absence of explicit solutions, which 
are not possible in general. We describe the generic evolution of 
the universe and show, in particular, that it is singularity free in 
the present class of models. Such models may stand on their own 
as interesting models for singularity free cosmology, and may be 
studied accordingly. They may also arise from string theory. We 
discuss critically a few such possibilities. 
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1. Introduction 



General Relativity (GR) is a beautiful theory and, more importantly, 
has been consistently successful in describing the observed universe. But, 
GR leads inevitably to singularities. In particular, our universe according 
to GR starts from a big bang singularity. We do not, however, understand 
the physics near the singularity nor know its resolution. The inevitability 
and the lack of understanding of the singularities point to a lacuna in our 
fundamental understanding of gravity itself. By the same token, a successful 
resolution of the singularities is likely to provide deep insights into gravity. 

In order to resolve the singularities, it is necessary to go beyond GR, 
perhaps to a quantum theory of gravity. A leading candidate for such a 
theory is string theory. However, despite enormous progress |T| , string theory 
has not yet resolved the singularities 0. 

Perhaps, one may have to await further progress in string theory, or a 
quantum theory of gravity, to be able to resolve singularities. However, there 
is still an avenue left open - resolve the singularities, if possible, within the 
context of a generalised Brans-Dicke model, and then inquire whether such 
a model can arise from string theory, or a quantum theory of gravity. 

Generalised Brans-Dicke model is special for more than one reason. It 
appears naturally in supergravity, Kaluza-Klein theories, and in all the known 
effective string actions. It is perhaps the most natural extension of GR [[|, 
which may explain its ubiquitous appearance in fundamental theories: Dicke 
had discovered long ago the lack of observational evidence for the principle 
of strong equivalence in GR. Together with Brans, he then constructed the 
Brans-Dicke theory |4[] which obeys all the principles of GR except that of 
strong equivalence. This theory, generalised in ||, has been applied in many 
cosmo logical and astrophysical contexts [fl|-[|IO|. 

Brans-Dicke theory contains a graviton, a scalar (dilaton) coupled non 
minimally to the graviton, and a constant parameter 10. GR is obtained 
when u = oo. In the generalised Brans-Dicke theory, referred to as graviton- 
dilaton theory, u is an arbitrary function of the dilaton || [J. Hence, it 
includes an infinite number of models, one for every function u. Thus, very 
likely, it also includes all effective actions that may arise from string theory 
or a quantum theory of gravity. 

This perspective logically suggests an avenue in resolving the singularities: 
to find, if possible, the class of graviton-dilaton models in which the evolution 
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of universe is singularity free, and then study whether such models can arise 
from string theory or a quantum theory of gravity. On the other hand, even if 
not deriveable from string theory or a quantum theory of gravity, such models 
may stand on their own as interesting models for singularity free cosmology. 
In either case, one may study further implications of these models in other 
cosmological and astrophysical contexts. Such studies are fruitful and are 
likely to lead to novel phenomena, providing valuable insights. 

Accordingly, in this paper, we consider a class of graviton-dilaton models 
where the function uj satisfies certain constraints. These constraints were 
originally derived in |TI| in a different context, and some of their generic 



cosmological and astrophysical consequences were explored in |TT], [12). Here 



we analyse the evolution of universe in these models which, as we will show, 
turns out to be singularity free. The essential points of the analysis and the 
results have been given in a previous letter [T3I . 



We consider a homogeneous isotropic universe, such as our observed one, 
containing matter. The relevent equations of motion cannot, in general, be 



solved explicitly except in special cases |Ll], [T2|]. Hence, a different approach 
is needed for the general case which is valid for any matter and for any 
arbitrary function uj, and which enables one to analyse the evolution and 
obtain its generic features even in the absence of explicit solutions. 

We will present such an approach below. We first present a general anal- 
ysis of the evolution and then apply it to describe in detail the evolution of 
universe in the present model. We show, in particular, that the constraints 
on uj ensure that the evolution is singularity free. The universe evolves with 
no big bang or any other singularity and the time continues indefinitely into 
the past and the future. 

An important question to ask, from our perspective, is whether a function 
uj as required in the present model, can arise from a fundamental theory. We 
discuss critically a few such possibilities in string theory. On the other hand, 
the present model may stand on its own as an interesting model for singularity 
free cosmology and may be studied accordingly. We mention some issues for 
future study. 

This paper is organised as follows. In section 2, we present our model. In 
section 3, we write its equations of motion in various forms, so as to facilitate 
our analysis. In section 4, we present the general analysis of the evolution. 
In section 5, we illustrate our method by applying the results of section 4 
to describe the evolution of toy universes and show that their evolution is 
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singularity free in the present model. In section 6, we describe the evolution 
of our observed universe and show that its evolution is singularity free in the 
present model. In section 7, we discuss further generalisations of our model. 
In section 8, we give a brief summary, discuss critically a few possibilities 
of our model arising from string theory, and mention some issues for future 
study. In the Appendix, we derive the sufficiency conditions for the absence 
of singularities which are used in the paper. 



2. Graviton-Dilaton Model 



We consider the following graviton-dilaton action in 'Einstein frame': [] 

S = — !— / d'x^g (-R + i(V0) 2 ) + S M (M, e~^% v ) , (1) 

where Gn is the Newton's constant, is the dilaton and ip{4>) ls an arbitrary 
function that charcterises the theory. Sm is the action for matter fields, 
denoted collectively by M.. They couple to the metric g^ v minimally and to 
the dilaton <fi through the function ip(4>). This function cannot be gotten rid 
of by a redefinition of g^ v except when the matter action Sm is conformally 
invariant - in that case, by definition, Sm(M, &~^^g^u) = Sm{M-,9hv) - 
which is assumed to be not the case here. 
We can define a metric 

g^u = e'^g^ , (2) 
and write the action ([!]) in 'Dicke frame': 

V=se* (-R + \ {^l - l) (V0) 2 ) + S M {M,g^) , (3) 



IottGat 



where ip^ = If i> has a finite upper bound tp < ip m ax < °o, then the 
factor e^ max can be absorbed into GV and the range of x/j can be set to be 
x/j < 0. Also, we will work in the units where Gn = e ™ x . 



Both forms of the action in (|I|) and (^) are, however, equivalent [|14 



In 'Einstein frame', the matter fields feel, besides the gravitational force, 
the dilatonic force also which must be taken into account in obtaining the 



_ i ft 2 

In our notation, the signature of the metric is ( — h ++) and R^\ T — - — 



2 dx^dx T 
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physical quantities. Whereas, in 'Dicke frame', the matter fields feel only the 
gravitational force and, hence, the physical quantities are directly obtained 
from the metric g^ v . On the other hand, equations of motion are often easier 
to solve if g^ v is used. 

Defining a new dilaton field x = e^, the action (|3]) can be written as 

uj(x) ,„ x 2 X 



167tGtv J \ X J 

where ou(x) is now the arbitrary function that characterises the theory. This 
is the form of the action used more commonly and, hence, we will also use 
it in this paper in most of what follows. Note that \ > since \ = & 1 . 
Furthermore, if x has a finite upper bound x < Xmax < °o, then the factor 
Xmax can be absorbed into Gn and the range of x can be set to be < x < 1. 
Also, we will work in the units where Gn = x ^ ax ■ 

Let Q = 2lo + 3. Then, (j) and "0(0) are related to x an d ^(x) as follows: 

d_X 

X 



= ip^ 2 or, equivalently, = f — \Jtt(x) ■ (5) 



Note that Q must be positive. If one starts from and "0(0), then the first 
relation gives Q in terms of 0. To obtain Q in terms of x-, '0(0) m ust be 
inverted first to obtain in terms of ip and then i/j = ln% must be used. 
If one starts from x an d ^(x) and uses x = e^, then the second relation 
gives in terms of which must then be inverted to obtain -0 in terms of 0. 
However, in the following, no inversion is necessary since no explicit forms 
for fi(x) and ip(4>) will be used. 

In the model we study here, the function ^(0) is required to satisfy the 
following constraints: 

lim if} = !|L n < - 



d n ip 
d(j) n 
dip 

# 



.W 3 
finite V n > 1, — oo < < oo 

at = only 

i2n 



lim ^(0) = -k<p ln , n > 1 , (6) 
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where fio > and k > are constants and n is an integer. The function 
ip(4>) is otherwise arbitrary. 

The first two constraints are important and were originally obtained in a 



different context in [11], 0], where some of their generic consequences were 
also explored. These constraints imply a finite upper bound on ip, i.e. ip{4>) < 
V'max < oo. Hence, the factor e^ max can be absorbed into Gn and the range 
of ip can be set to be ip < 0. Thus there is atleast one critical point where 
# = 0, and V = 0. 

The third constraint is not important, but assumed here for the sake of 
simplicity only. It implies that there is only one critical point, a maximum, 
of if}((f)). As explained above, we can take VWs = 0. Thus, the range of ijj is 
given by — oo < ip < and, hence, that of x — e by < x ^ 1- Also, by 
adding a suitable constant to (p, we can take the maximum of i/j(4>) to occur 
at (f> = with no loss of generality. 

The fourth constraint follows from the Taylor expansion of i/j(4>) near its 
maximum at <ft = 0. Note, as follows from @, that Q —>■ oo as (ft 0. 

The constraints in (|6]), expressed equivalently in terms of \ and 
become 

O(0) = 

- — = finite V n > 1, < y < 1 

dx n 

— > oo at x = 1 only 

1 



hmfi = ^<a<l, (7) 



where f2 > and Qi > are constants, and the range of a follows naturally 
from (H) and (|6]). We further assume that Q(x) is a strictly increasing func- 
tion of x- This constraint is not important, but assumed here for the sake of 
simplicity only. The function fl(x) is otherwise arbitrary. 

Note, for later reference, that Q — > oo and = — fpr(l ~ x) 4 "" 1 i n the 
limit x ~ * 1- Hence, for the range of a given in (|7|), which follows naturally 
from (|6]), we have that - > in the limit ^ — > 1. This ensures, as 
will become clear later, that the present model satisfies the observational 
constraints imposed by solar system experiments. 

A vast class of functions equivalently exists satisfying the 

above constraints. A simple example is Q = Qq+Qi((1 — x)~ 2q ~ !)• However, 
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no explicit form for ft(x) wm be used in the following. 



3. Equations of Motion 

In the following, we take the "matter" to be a perfect fluid with density 
p and pressure p, related by p = 7/9 where — 1 < 7 < 1. Note that the 
value of 7 indicates the nature of "matter". For example, 7 = 0, |, or 1 
indicates that the "matter" is dust, radiation, or massless scalar field respec- 
tively. Hence, we enclose "matter" within inverted commas where necessary 
in order to remind the reader that "matter" does not mean dust only. In the 
following, we study the evolution of a flat homogeneous isotropic universe. 
(The following analysis can, however, be extended to curved universes also 
and the main results remain unchanged.) The line element is given by 

ds 2 = -dt 2 + e 2A V(dr 2 + r 2 dS 2 ) , (8) 

where e A is the scale factor and 0LS2 is the line element on an unit sphere, 
and the fields depend on the time coordinate t only. The equation of motion 
for A, following from (f|), is given by 

x x 2 Qx 

where upper dots denote t-derivatives. Solving equation @ for A, we obtain 



2 X V QX 12x 2 



where e = ±1 and we have used f2 = 2u + 3. The square roots in (|iy) and 
in the following are to be taken with a positive sign always. The remaining 
equations of motion, following from ([|), are given by 

X + 3ix+^ = (1 " 37)P (11) 

P = Poe- 3(1 ^ )A , (12) 

where we have used p = jp and po ^ is a constant. Equations of motion 
following from (ffl) also contain the equation for A: 



A = 2A + 3A 2 + ^ + ^ + * + ^ = 0. (13) 
X %X X 2x 
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In general, this equation follows from (^), (|TT|), and ( |T2"D also. But there is 
a subtlety pointed out in What one obtains from (|), ([ll]), and (|l^) is 
the following: [] 

( i+ |^M = ' (14) 

where the expression for A is defined in (|13D. Equation ^ |X3|) then follows if 
and only if (A + 7^ 0, i.e. if and only if e A is not proportional to \ 
This subtlety should be kept in mind. However, as checked explicitly, the 
solutions presented in this paper all satisfy equations @ - (|13|). 

Equations (pi]) and ([ll]) can also be written in different forms. Equation 
( |TTD can be integrated once to obtain 

— 3 A 

X(t) = °(<,(t)+c) (15) 



where „ {t) = t 1 - 3 ^ f' dt ^ (16) 



and c = xv'fie 3 ' 4 |t=t; is a constant. An equation relating A and x can now 
be derived. Divide equation (|T0| ) by ^ and substitute ( |T5| ) for We then 
have 

<*>^. <»> 

All the above equations can be equivalently written in terms of <fi and 
■)/>(</>) also. For example, equations (|T0|) and ( |TTD become 



A+- = el + — (18) 

2 V 6 12 v ; 

+ 3A0 + ^> = (1 ~ 2 37)P ^ , (19) 



2 Although the derivation of ( J14| ) is straightforward, an outline of the steps involved 
is perhaps helpful. In equation (flf), restore the p-dependence by replacing (1 — 37)p by 
(p — 3p). Now, differentiate (^) and use p = — 3A(p + p). Then substitute (^) for and 

( pi] ) for ci), equivalently f2. Again substitute (j9j) for . The result is equation ( |l4| ) p"5[ ]. 
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while equation (|T5| ) and ( |T6D become 



where 



4>(t) 

a(t) 



-3A-tp 



>(t) + C) 



;i - 37)po 



dt t/j^e 



-37 A 



(20) 
(21) 



and c = 0e 3A+ ^ | t=ti is a constant. 

We now make a few remarks. Firstly, under time reversal t — > —t, e — > — e 
in equations ([10|) and (|18D, whereas equations ( |TTD and (|I9| ) are unchanged. 
Thus, evolution for — e is same as that for e, but with the direction of time 
reversed. 

Secondly, in general, there will appear in the solutions positive non zero 
arbitrary constant factors e A °, po, and xo m front of e A , p, and x respectively. 
However, it follows from (|ToD and ([□]) that they can all be set equal to 1 with 

no loss of generality if time t is measured in units of y l xoe 3 ^^ o_ ^ r pj iere f or g ) 

in the solutions below we often assume that t is measured in appropriate 
units and, hence, set these constants equal to 1. Note, however, that the 
model dependent constants associated with fl(x), such as Q or Qi in (0), 
cannot be set equal to 1. 

Thirdly, if p = or 7 = | then a(t) = in (Jl5|). If it is also possible 
to perform the integrations given below then solutions to equations fll0|) and 
( ITS] ) can, in principle, be obtained in a closed form. To see this, define two 
new variables t and x as follows: 



dt = y/xdt 



(22) 



Note that these variables appear naturally in 'Einstein frame' given in ([I]) 
and (0). After some algebra, equations (|Il]) and ( [T3|) become 



dA 
~dl 

ldx 
X di 

from which it follows that 

xdA 
dx 



ece 



-3A 



ce 



12 

3A 



1 + 



2p e 



2.1 



n / + 2p^ 



\12 



(23) 
(24) 

(25) 
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In the above equations e = ±1 as in ( |T0"D and c is a constant as in flI5|) . 

It turns out that A-integrations in equations ( p3|) and ( p5|) can be per- 
formed explicitly. Omitting the details, which are straightforward, the result 
is as follows. When p and, hence, p = 0: 

When p ^ 0, define = x fW e A . Then 



Wl + y 2_ lm y 1 + 2/ 2 + ?/ ) = "75 Hr ( 28 ) 




lny-huVl + ^ + l) = ^ (29) 



One thus has an explicit relation between A and t, and also between A 
and x if one can obtain / ^*\/VL in a closed form. Thus, one has x(t). If 



one can obtain / —j= also in a closed form then one has an explicit relation 

between t and t. Therefore, in principle, one has A(t) and x(i) in a closed 
form. 

Evidently, solutions can be obtained thus in a closed form only in special 
cases. Even then the details of the required integrations and inversions of 
functions are likely to obscure the general features of the solutions. Such 
explicit solutions have been studied in [TTJ for a few specific functions fl{x), 
and in |12| for any arbitrary function fl(x) but with p = 0. In the present 



paper, however, we follow a different approach which is valid for any matter 
and for any arbitrary function O(x)- 

4. Analysis of the Evolution 



Our main goal in this paper is to determine whether the evolution of the 
universe in the present model is singular or not. The task is trivial if one 
can solve the equations of motion (|T^) - flT2] ) for arbitrary functions £l(x) ( or 
il>((j))). However, explicit solutions can be obtained only in special cases as 
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described above. Hence, a different approach is needed for the general case 
which is valid for any matter and for any arbitrary function fl(x), and which 
enables one to analyse the evolution and obtain its generic features even in 
the absence of explicit solutions. 

We will present such an approach below. We analyse the evolution and 
obtain its generic features for the general case where matter and f)(x) ( or 
ip(4>)) are arbitrary. We will show, in particular, that when O(x) (or ifj(4>)) 
satisfies the constraints in (^) (or (|6])), the evolution of the universe, such as 
our observed one, is singularity free at all times. 

To prove the presence of a singularity it suffices to show that any one of 
the curvature invariants, usually the curvature scalar, diverges. However, to 
prove the absence of a singularity, it is necessary to show that all curvature 
invariants are finite. As proved in the Appendix, a sufficient condition for 
all curvature invariants to be finite is that the quantities in fl5"3"| ) be all finite. 
Some or all of these quantities have a potential divergence when, for example, 
e A and/or x vanish or when Q — ► oo. Hence, in the following analysis, we 
pay particular attention to these quantities, determining their finiteness or 
otherwise. We perform our analysis, as far as possible, in terms of x and 
Q(x). In the limit as x ~ > 1> however, it is sometimes necessary to perform 
the analysis in terms of <fi and ip{4>), as will become clear below. 

To proceed, we need initial conditions at an initial time, say t = U. That 
is, we need A(ti), x(U), and x(U) or equivalently Q(ti). Here and in the 
following Q(t) means Q(x{t)), i.e. the function fl(x) evaluated at time t. 
Note that if x(ti) — and x(ti) = 1> then Q(ti) = oo and the evolution 
is same as in Einstein's theory. In particular, there is a singularity at a 
finite time in the past. Therefore, we assume that the initial conditions in 
the following refer always to generic values of A(ti), x(U), and fl(tj) only, 
i.e. A{ti) and are finite and non infinitesimal and < oo in the 

following. 

The analysis will be carried out for any arbitrary function fi(x); satis- 
fying only the constraints in (|7]). Hence, only the generic features of the 
evolution can be obtained but not the numerical details which require choos- 
ing a specific function for fi(x)- Therefore, precise numerical values for A(U) 
and x{ti) are not needed. Their signs alone are sufficient. We also need to 
know whether or not fi(£j) ^ 3, and also the value of e in fliCf). Thus, there 
are 16 sets of possible initial conditions: 2 each for the signs of A(ti), x(U), 
and e, and 2 for whether Q(ti) £ 3 or not. 
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However, it is not necessary to analyse all of these 16 sets. The evolution 
for — e is same as that for e, but with the direction of time reversed. Hence, 
only 8 sets need to be analysed. Now, let e = 1. If x < then A can not 
be negative for any value of Q, see fllPf) . Similarly, if x > and Q > 3 then 

also A cannot be negative since \J -§^ + > ^ > in equation (|10D . 
Thus, there remain only 5 sets of possible initial conditions. We choose these 
sets so as to be of direct relevence to the evolution of observed universe and 
analyse each of them for t > ti. In the course of the analysis, equivalent 
forms of the equations of motion given in section 3 will become useful. We 
then use these results in sections 5 and 6 to describe the generic evolution of 
the universe. 

Two remarks are now in order. First, the amount and the nature of 
the dominant "matter" in the universe are given by po an d 7 where — 1 < 
7 < 1. The later varies as the scale factor e A of the universe evolves. If e A is 
increasing, i.e. if A > 0, then the universe is expanding, eventually becoming 
dominated by "matter" with 7 < ~, when such "matter" is present. For the 
observed universe, which is known to contain dust for which 7 = 0, one may 
take 7 < 0. Thus, by choosing ti suitably we may assume, with no loss of 
generality, that 7 < | (< for observed universe) if A > 0. Similarly, if 
A < then we may assume, with no loss of generality, that 7 > | when such 
"matter" is present. This is valid for the observed universe which is known 
to contain radiation for which 7 = |. 

Second, consider a{t) in (|T^) or (|2"0"|). The integrand is positive and oc 
(The dependence on p can be absorbed into the unit of time, as explained in 
section 3.) Hence, the value of the integral is controlled by Q. For example, 
in the limit \ — > 1, it is controlled by the constant Qi in (^|). We will use this 
fact in the following. Also, since t > t^, the sign of a(t) is same as that of 
(1 — 37). Moreover, a{t) = if po = or 7 = | corresponding to an universe 
containing, respectively, no matter or radiation only. 

4 a. A{ti) > 0, x(U) > 0, and fi(t«) > 3 



These conditions imply that e = +1 in ( TOf) and the constant c > in 
(|l5|). For t > U, the scale factor e A increases and, eventually, 7 can be taken 
to be < I when such "matter" is present. In fact, 7 < for the obeserved 
universe. Then, (1 — 37) > and, hence, a(t) increases. Thus, for t > U, 
X > and x an d ^ both increase. Since Q > 3, it follows from equation 
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( |IOD that A > for t > U and, hence, e A increases. Eventually, x ~ > 1 an d 

q ^ n x (i - x)- 2 «. 

In the limit as x — > 1, (<r(£) + c) > is a constant to an excellent 
approximation. Then, equation (O) can be solved relating x an d e A : 



e A = (constant) (1 - x) • (30) 



Substituting this result in equation (|T5|) then yields the unique solution in 
the limit x ~ > 1 ; 



,A _ „An 



e A = e^t 3 ^) , x = 1 - Xot ( - 1 +^<- 1 -^ , (31) 
where A and Xo > are constants and t is measured in appropriate units. 



Since ^ — ► 1, it follows from ([H|) that t — > oo. It can now be seen that 
(cr(t) + c) > is indeed a constant to an excellent approximation. Also, the 
quantities in (|53|) are all finite, implying that all the curvature invariants are 
finite. Thus, there is no singularity as x - ¥ 1 and, thus, for t > U. 

Note that the above relation between x and e A and, hence, the solution 
( pip is not valid if p = or 7 = 1 i.e. if the universe contains no matter 
or "matter" with 7=1 only. Although these cases are not relevent for the 
observed universe, we will now consider them for the sake of completeness. 
The results, however, are useful in describing the evolution for toy universes 
in section 5. 

p = 

p = implies that po = 0. In this case, a(t) = identically and the 
analysis was first given in ]T2|, albeit differently. The dynamics turns out to 



be clear in terms of <f> and ifr((f))- These variables, although equivalent to x 
and fi(x); are particularly well suited in describing the evolution near = 0, 
or equivalently x = 1, where Q — > 00. In fact, in this limit, it is not possible 
sometimes to perform the analysis in terms of x and f2(x). 

At £j, we have A > 0, <p > and, with our normalisation, <p < 0, r/> < 
and ipf > 0. Also, t/>(0) has a maximum -i/Wz = at = 0. The equations 
of motion flI5|) and flT9|) are now given by 



A = (-V3VV + sign(0)) ^= , = ce~ 3A ^ , (32) 
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where c > 0. For t > U, the scale factor e A and continue to increase. 
Eventually, — > 0. In this limit, -0 — > and ^ — > 0. The solution to ( j32l) is 
then given by 

e A = e A H \ ^ = c l n £ , (33) 

where £ is measured in appropriate units. From these expressions, it is clear 
that the field crosses the value at t — 1 (in appropriate units) and 
becomes positive, while the scale factor e A remains finite and increasing. For 
> 0, ip((f) decreases and ip<i> < 0. Equivalently, x an d &(x) both decrease. 
Thus, we now have A > 0, x < 0, and e = 1. Further evolution is analysed 
in section 4b. 

Note the crossover at = where crosses and becomes positive. In 
terms of x an d ^(x), this amounts to the following: x reaches 1, correspond- 
ingly Q reaches oo, and then x "reflects" back and decreases, correspondingly 
Q decreases. This crossover /reflection takes place in a finite time and with 
e A remaining finite during the process. It is this phenomenon of crossover 
that can be seen only in terms of and ip{4>), and not x an d ^(x)- 

7 = 1 

In this case, the universe contains "matter" with 7 = 1 only. Hence, p = 
Poe~ 6A and a(t) is negative. However, as x ~~ * 1) a (t) oc ^7=, see equation 
(0), and remains finite as can be verified a posteriori. Also, the magnitude 
of a(t) can be made as small as necessary by choosing the model dependent 
constant Qi sufficiently large. Hence, as x — > 1 5 the factor (a(t) + c) can be 
taken to remain positive and constant. 

Now, VLx 2 = (constant) e~ 6A and its dependence on e A is same as that of 
p. Therefore, in the limit x ^ 1 or equivalently Q — > 00, the solution to (|i~8|) 
and (|20|) and, hence, the evolution turn out to be identical to the p = case. 
In particular, the dynamics is clear in terms of and only. The field 
crosses the value (at t = 1 in appropriate units) and becomes positive, 
while the scale factor e A remains finite and increasing. This crossover cannot 
be seen in terms of x an d ^(x)- 

For > 0, '0(0) decreases and -0^ < 0. Equivalently, x an d ^(x) both 
decrease. Thus, we now have A > 0, x < 0, and e — 1. Further evolution is 
analysed in section 4b. 

4 b. A(U) > 0, x(U) < 0, fl(tj) > 3, and e = 1 
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The constant c < in (|1"5|). For t > ti, x an d ^ decrease and the scale 
factor e A increases. Equation now becomes 

dx 

where K{t) is given in (p!7|). 

Consider first the case p = or 7 = 1, discussed in section 4a and which 
led to the present phase. Then, p (l — 37) < and, therefore, a(t) < 0. 
Thus, + c) < c < and, hence, x(t) is negative and non zero, implying 
that decreases for t > t^. Also, K{t) remains finite and it follows from 
equation ([17|) that ^ < 0. Thus, ^4 > since x < and, hence, e A increases 
for t > ti. 

Eventually, as t increases, x ~~ > 

and e A -> 00. Also, 2 ^ ( g+g < 1 m 
( p7|) since po = or 7 = 1. Equation (|T7| ) can then be solved relating x an d 



e A : 



e = (constant) x 1 > (34) 

where x — >■ and f2 is given in (0). Substituting this result in equation (|15|) 
then yields the unique solution in the limit x ^ 0: 

e A = e A °(t-t o r , X = Xo(t-t ) m , (35) 

where A , Xo, an d to > U are some constants and 



s + ym~ -2 

n = — ; , m= ; . 36 

3(1 + V^h) 1 + V^h 



Note that m < 0. Then, since x ~ ^ 0, it follows from flHoD that t — > 00. 

Thus, as x — * 0, £ — > 00 and e A — > 00. It can also be seen that the 
quantities in (|53f) are all finite for ti < t < 00, implying that all the curvature 
invariants are finite. Thus, there is no singularity for U <t < 00. 

Consider now the case where p and 7 are non zero and have generic values. 
We have Afa) > and x{U) < 0. Hence, e A increases and x decreases. As 
e A increases, 7 can eventually be taken to be < |, when such "matter" is 
present. In fact, 7 < for the obeserved universe. In such cases where 7 < 0, 
a(t) grows faster than t as can be seen from fll6|) . Then, as t increases, the 
factor (er(t) + c) becomes positive. 
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This implies that \i initially negative, passes through zero and becomes 
positive. We then have A > and x > 0. As t increases further, x and, 
hence, £l(x) increase. If Q > 3 then further evolution proceeds as described in 
section 4a. If Q < 3 then, upon making a time reversal, the initial conditions 
become identical to the one described in section 4e. The evolution in the 
present case then proceeds as in section 4d for the case which leads to the 
initial conditions of section 4e. 

4 c. A(ti) < 0, x(ti) < 0, and Slfc) > 3 

These conditions imply that e = — 1 in fllOl ) and the constant c < in 
(|i~5|). For t > U, the scale factor e A decreases and, eventually, 7 can be taken 
to be > I when such "matter" is present. In fact, this is the case for the 
observed universe. Then (1 — 37) < and a(t) decreases or remains constant. 
Therefore, the factor (cr(t) +c) < c < and, since e A decreases, we have that 
x(t) < x{ti) < for t > t{. Thus x and, hence, fl decrease for t > t; L . Also, 
^ > since x < and A < 0. 

In (|i~7D , K{ti) > 1 since Q > 3. From the behaviour of e A , %, and 
(cr(t) + c) described above, it follows that K decreases monotonically. The 
lowest value of K is ^ < |, achievable when x vanishes. This, together with 
K(ti) > 1 and the monotonic behaviour of K(t), then implies that there 
exists a time, say t = t m > U where K{t m ) = 1 with x(t m ) > 0. Hence, 
^(tm) = 0. Therefore, A(t m ) = since x(t m ) is non zero. This is a critical 
point of e A and is a minimum. Also, equation ([[7]) can be written as 

A(ti) ~ A(t m ) = f XiU) p (-1 + VK)= finite , (37) 

J X(tm) ^X 

where the last equality follows because both the integrand and the interval 
of integration are finite. Therefore, A(t m ) is finite and > — 00 and, hence, 
e A{t m ) j g g n jt e anc [ non vanishing. (] 

Thus, for t > ti, the scale factor e A continues to decrease and reaches a 
non zero minimum at t = t m . The precise values of t m , A(t m ), x(t m ), and 
fl(t m ) are model dependent. Hence, nothing further can be said about them 

3 The existence of such a non vanishing minimum of e A has also been deduced in the 
variable mass theories of Q. However, singularities may subsequently arise in these theo- 
ries. See section 7 also. 
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except, as follows from K(t m ) = 1, that £l(t m ) £ 3 in general and = 3 when 
Po = 0. 

However, the above information suffices for our purposes. It can now be 
seen that the quantities in (^) are all finite, implying that all the curvature 
invariants are finite. Thus, there is no singularity for ti < t < t m . 

For t > t m , one has A(t) > 0, x(t) < 0, and fl(x(t m )) < 3 by continuity. 
Further evolution is analysed below. 

4 d. A(ti) > 0, x(ti) < 0, < 3, and e = -1 

This implies that the constant c < in (|15D and that < in (|T?[). 

Hence, K(ti) < 1. For t > U, the scale factor e A increases and, eventually, 7 
can be taken to be < | when such "matter" is present. In fact, 7 < for the 
observed universe. Then, (1 — 37) > and, hence, a(t) increases. Thus, e A 
is increases, x decreases, and (<r(t) + c) increases. Now, depending on p , 7, 
and the details of K(t) in (|l~7|) may or may not remain < 1 for t > ti. 

Consider first the case where K(t) remains < 1 for t > ti. It then follows 
that (<j(t) + c) cannot have a zero for t > ti and, since c < 0, must remain 
negative and non infinitesimal. Note that 7 must be > for otherwise a(t) 
grows as fast as or faster than t, making (a(t) + c) vanish at some finite time. 
Let 

lim (a(t) + c) = c e (38) 

t — >oo 

where c e is a negative, non infinitesimal constant. We then have for <t < 
00, c < (cr(t) + c) < c e < and, hence, x{t) < 0- Therefore, eventually 
X{t) -> 0. Since if(t) < 1 and x{t) < 0, it also follows from (0) that ^ < 
and, hence, A{t) > 0. We will now consider the limit x ~~ > 0. 
K(t) < 1 implies that as^— > 0, fC^^<l where 

fie^o + ^ limxe 3 ^^ (39) 

is a constant and Qo is given in (^). Now, as x ~ ¥ 0, equation (|17D can be 
solved relating x an d eA '- 

e = (constant) x 6 • (40) 
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Note that (3 — > since Q e < 3. Hence, e A — > oo as x 0. 



Substituting ( fiO|) in equation (fL5l) then yields the unique solution in the 
limit x — > 0: For fi e 7^ ~, 



e = e (t - sign(m)t) , x = Xo (<o - sign(m)t) , (41) 

where t is measured in appropriate units, A , Xo> an d > ^ are some 
constants, and 



n 



y/3TT e -2 



3(1 - V3H^ ' 1 - v 7 ^ 



m = r = . (42) 



For n e = §, 



where c e is defined in 

Let Vt e > \. Hence, m > 0. Also, n < because f2 e < 3. Since x ~ y 0, it 



e A = e Ao e -^ f x = Xoe Ce(t-t ) f (43) 



follows from (|4l| ) that t ^ to > U, which implies that x vanishes at a finite 
time to- in this limit, the scale factor e A — > 00 since n < 0. The quantities 
in ( j53|) ) for example 2£, also diverge, implying that the curvature invariants, 
including the Ricci scalar, diverge. Thus, for D, e > |, there is a singularity 
at a finite time to- 

Let f2 e < |, which is the case in our model, see (0). When Q e < |, m < 0. 
Also, n > because f2 e < 3. Since x ~ > 0, it follows from (|4l| ) that t — > 00. 



In this limit, the scale factor e A — > 00 since n > now. The same result 
holds also for the solution in (f43"D with tt e = ~. 

For -KT(t) to be < 1, (cr(t) + c) must not be too small. Since c < 
and a(t) > 0, it necessarily implies that a(t) must remain finite. It follows 
from the above solutions that a(t) can remain finite only if -y/3fi e > -jz^f-, 

equivalently, 7 > Under this condition, it can be checked easily that 

lim x ^ xe 3 ^ 1-7 ^ = 0. Hence, f2 e = fi as follows from ([39D . 

It can now be seen, for Q e = Q < |, that the quantities in ( p3| ) are all 
finite for U < t < 00, implying that all the curvature invariants are finite. 
Thus, there is no singularity for t^ < t < 00 when Q < §■ 

Consider the second case where K(i) in (fl7|) may not remain < 1 for all 
t > ti. This is the case for the observed universe where 7 < and, hence, 
a{t) grows as fast as or faster than t. Then there is a time, say t — t%, when 
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the value of K = 1. It follows that c < (cr(ti) + c) < and, hence, x(t) < 
for U < t < t\. Also, x(ti) is non vanishing since the case of x(ti) — > is 
same as the case described above. Therefore, we have that ^p(ii) = which 

implies, since x(*i) 7^ 0> that ^(^i) = 0- This is a critical point of e A and is 
a maximum. 

For t > t\,K(t) becomes > 1 and, hence, A(t) < 0. Also, x < and 
Q(t) < 3. Then e = — 1 necessarily. Further evolution is analysed below. 

4 e. A(ti) < 0, < °> and O(ti) < 3 

The above conditions imply that e = — 1, the constant c < in 
jfc{ti) > 0. Hence, K(ti) > 1. For t > ti, (cr(t) + c) which is negative at i; 
may or may not vanish for non zero X- 

Consider the case where (cr(t) + c) does not vanish and remains negative 
for non zero x- Therefore, x(£) < and, hence, x(t) decreases. For t > U, 
the scale factor e A decreases and, eventually, 7 can be taken to be > ~ when 
such "matter" is present. In fact, this is the case for the observed universe. 
Then, (1 — 37) < and a(t) decreases. Hence, (cr(t) + c) also decreases. 

It follows from ( |TTD that K, initially > 1, is now decreasing since x an d 
e A are decreasing and (cr(t) + c) 2 is increasing. Its lowest value = ^ < |, 
achieveable at x = 0. Therefore, there exists a time, say t = tM, where 
K = 1 and, hence, ^ = 0. Clearly, x(^m) and x(£m) are non zero for the 

same reasons as in section 4c, implying that A(tjy) = 0. This is a critical 
point of e A and is a minimum. 

The existence of this critical point of e A can be seen in another way also. 
As (cr(t) + c) < grows in magnitude, it follows from equation flTJD that 
f^x 2 oc e _6A , whereas p is given by (p!2|). Note that x(t) < and, hence, x(t) 
decreases and eventually — > 0. Then, in the limit x ~^ 0> the term in 
equation ( p!Q[) dominates ^ term. Hence, in this limit where — > S7 < ^ < 3, 

we have that A > 0. Since A < initially, this implies the existence of a zero 
of A. This is a critical point of e A and is a minimum. 

Note that all quantities remain finite for U < t < tj^- In particular, the 
quantities in ( |53D are all finite, implying that all the curvature invariants are 
finite. Thus, there is no singularity for U < t < tj^- 

For t > t^j, we have A > 0, x < 0, f2 < 3, and e = —1. Further 
evolution then proceeds as described in section 4d. The evolution can thus 
become repetetive and oscillatory. 
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Consider now the case where (cr(t) + c) vanishes, say at t — t m /, for 
non zero X- Hence, x(t m ') vanishes and this is a minimum of X- F° r t > 
t m r, (cr(t) + c) and x are positive by continuity and, hence, x increases. Also, 
A < and the scale factor e A decreases. Eventually, 7 can be taken to be 
> I when such "matter" is present. In fact, this is the case for the observed 
universe. Then, (1 — 37) < and a(t) decreases. Hence, (cr(t) + c) also 
decreases. 

Now, as x increases and — > 1, (a(t) + c) may or may not vanish with 
X < 1- If (&(t) + c) vanishes with % < 1, then x vanishes. This is a critical 
point of x an d is a maximum, beyond which we have A < 0, x < 0, and 
e = — 1. Further evolution then proceeds as described in section 4c if Q > 3, 
and as in section 4e if f2 < 3. The evolution can thus become repetetive and 
oscillatory. 

If (a(t)+c) does not vanish and remain positive for x < 1 then, eventually, 
X — > 1. It follows from equation ([15]) that fix 2 °^ e~ 6A , whereas p is given 
by (12). Note that A < and e A is decreasing. Then, in this limit, the 
term in equation ( pf) dominates or, if 7 = 1, is of the same order as the £ 
term. 

Thus we have A < 0, x > 0, x — *" 1 ; equivalently — > 00, and e = — 1. 
Also, term in equation ( Jl~0| ) dominates or, if 7 = 1, is of the same order 
as the ^ term. But this is precisely the time reversed version of the evolution 
analysed in section 4b and in section 4a for the case which led to the initial 
conditions of section 4a, where the dynamics is clear in terms of <fi and i/){(f))- 
Applying the results of sections 4a and 4b, it follows that <fi will cross the 
value after which and, hence, \ begins to decrease. We then have A < 0, 
X < 0, and Q > 3. Also, e = — 1 necessarily. Further evolution then proceeds 
as described in section 4c. The evolution can thus become repetetive and 
oscillatory. 



5. Evolution of toy Universes 



We now use the results of the analysis in section 4 to describe the generic 
evolution of universe in the present model. For the purpose of illustration we 
first describe in this section the generic evolution of three toy universes: One, 
where the universe contains no matter, i.e. p = 0. Two, where the universe 
contains radiation only, i.e. p 7^ but 7 = |. And three, where the universe 
contains massless scalar field only, i.e. p 7^ but 7=1. Clearly, however, 
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none of these toy models are realistic. The generic evolution of a realistic 
universe, such as our observed one, will be described in the next section. 

Note that a(t) = in the first two models. As described in section 3, it 
may then be possible to obtain closed form solutions, although only in special 
cases fi"0| , Moreover, the details of the specific calculations obscure the 
general features of the solutions. Hence, and also for the purpose of illustra- 
tion, we describe the evolution of the toy universes within the framework of 
our analysis. However, the explicit evolutions described in |L0|, |12[ for specific 
cases all conform to the ones described by the present general analysis. 

We first describe the evolution for t > U taking, as initial conditions 
A{ti) > 0, x(U) > 0, and Q(ti) > 3. Then, e = 1 necessarily. To describe 
the evolution for t < ti, we reverse the direction of time and take, as initial 
conditions, A(ti) < 0, x{U) < 0, and fi(tj) > 3. Then, e = — 1 necessarily. 
The required evolution is that for t > ti in terms of the reversed time variable, 
which is also denoted as t. The results of section 4 can then be applied 
directly. 

When the initial conditions are different, the evolution can again be anal- 
ysed along similar lines as follows. However, the main result that the evolu- 
tion is singularity free remains unchanged. 

5 a. p = 0, t > ti 



The evolution for this case was first described in JL2f where explicit so- 
lutions to the equations of motion were obtained for any arbitrary function 
in terms of the variables defined in (|22|) . We will now describe this 
evolution within the framework of the analysis presented in section 4 which 
is valid for any matter and for any arbitrary function fi(x)- 

Initially, we have A(ti) > 0, xfti) > 0> an d Q{ti) > 3. Then, e = 1 
necessarily. For t > ti, the evolution proceeds as described in section 4a. 
Both e A and \ increase. Hence, Q increases. Eventually, x ~ > 1 an d ^ — ► oo. 
The dynamics near x = 1, equivalently <p = in our normalisation, is clear in 
terms of <fi and ip ((/)). It follows from the analysis of section 4a that cf) crosses 
and continues to increase. Then, ip{4>) and, hence, Q decrease. 

We then have A > 0, but x < 0. The evolution then proceeds as described 
in section 4b. e A increases whereas x an d, hence, fl decrease. As t — ► oo, e A 



and x evolve as given in ([35|) and (j3~6P . 
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In particular, it can be seen that the quantities in (|53| ) are all finite for 
ti < t < oo, implying that all the curvature invariants are finite. Hence, the 
evolution is singularity free. 

t < U, equivalent ly (— t) > (-U) 

To describe the evolution for t < ti, we reverse the direction of time. 
Then, in terms of the reversed time variable also denoted as t, we have 
A(ti) < 0, x(ti) < 0, and > 3. Then, e = — 1 necessarily. For t > ti, 

the evolution proceeds as described in section 4c. Both e A and x decrease. 
Hence, Q decreases. Then, as shown in section 4c, there exists a time, say 
t = t m > U where K{t m ) = 1 in equation ([T7|), e A( * m) > 0, x(£ m ) > 0, and 
x(tm) < 0. Note that, in this case where p and, hence, po vanish, we have 

Q(tm) = 3. 

As shown in section 4c, K(t m ) = 1 implies that A(t m ) = 0. Hence, the 
scale factor e A reaches a minimum. For t > t m , we then have A > 0, x < 0) 
Q < 3, and e = — 1. The evolution proceeds as described in section 4d. 

It follows that a(t) = since p and, hence, po vanishes. Hence, in equation 
(|38|), c e = c which is negative and non infinitesimal. Therefore, as described 
in section 4d, e A increases whereas x an d, hence, Q decrease. As t — » oo, e A 
and x evolve as given in ( P]) - (^3|). 

In particular, it can be seen that, for Q,q < |, the quantities in (^) are 
all finite for U <t < oo, implying that all the curvature invariants are finite. 
Hence, the evolution is singularity free. 

Thus, the evolution of the toy universe with p = and with the initial 
conditions A(ti) > 0, x{U) > 0, and Q(ti) > 3 proceeds in the present model 
as follows. For t > t{, the scale factor e A increases continuously to oo. The 
field x increases, reaches the value 1, and then decreases continuously to 0. 
Correspondingly, Q increases, reaches oo, and then decreases continuously to 
tt , given in (0). 

For t < ti, the scale factor e A decreases, reaches a non zero minimum, 
and then increases continuously to oo. The field x decreases continuously to 
0. Correspondingly, fl decreases continuously to Q , given in (0). 

Also, the curvature invariants all remain finite for — oo < t < oo. Hence, 
in the present model, the evolution of the toy universe with p = is singu- 
larity free. 
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When the initial conditions are different, the evolution can again be anal- 
ysed along similar lines. However, the main result that the evolution is 
singularity free remains unchanged. 

5 b. 7 = |, t > U 



The evolution for this case was first described in [ 10) where explicit so- 



lutions to the equations of motion were obtained for a few specific functions 
£l(x)- We will now describe this evolution within the framework of the anal- 
ysis presented in section 4 which is valid for any matter and for any arbitrary 
function Q(x)- 

Initially, we have A(ti) > 0, xiU) > 0, and Q(ti) > 3. Then, e = 1 
necessarily. For t > t iy the evolution proceeds as described in section 4a. 
Both e A and x increase. Hence, Q increases. As t — > oo, e A and x evolve as 
given in ( |3~TD where 7 = | now. 

In particular, it can be seen that the quantities in (|53|) are all finite for 
ti <t < oo, implying that all the curvature invariants are finite. Hence, the 
evolution is singularity free. 

t < U, equivalently (—t) > (-U) 

To describe the evolution for t < ti, we reverse the direction of time. 
Then, in terms of the reversed time variable also denoted as t, we have 
A(ti) < 0, x(ti) < 0, and Q(ti) > 3. Then, e = — 1 necessarily. For t > ti, 
the evolution proceeds as described in section 4c. Both e A and x decrease. 
Hence, Q decreases. Then, as shown in section 4c, there exists a time, say 
t = t m > U where K{t m ) = 1 in equation (0), e A{tm) > 0, x(t m ) > 0, and 
X{tm) < 0. Note that tt(t m ) < 3. 

As shown in section 4c, K(t m ) = 1 implies that A(t m ) = 0. Hence, the 
scale factor e A reaches a minimum. For t > t m , we then have A > 0, x < 0> 
Q < 3, and e = — 1. The evolution proceeds as described in section 4d. 

It follows that a(t) = since 7=1- Hence, in equation fl38|) , c e = c 
which is negative and non infinitesimal. Therefore, as described in section 
4d, e A increases whereas x an d, hence, Q decrease. As t — > oo, e A and x 
evolve as given in ([y]) - d43f) . 

In particular, it can be seen that, for < ~, the quantities in (^) are 
all finite for U <t < oo, implying that all the curvature invariants are finite. 
Hence, the evolution is singularity free. 
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Thus, the evolution of the toy universe containing only radiation, i.e. 
p 7^ but 7=5, and with the initial conditions A(ti) > 0, x(U) > 0, and 
Q(ti) > 3 proceeds in the present model as follows. For t > t i: the scale 
factor e A increases continuously to 00. The field x increases continuously to 
1. Correspondingly, fl increases continuously to 00. 

For t < ti, the scale factor e A decreases, reaches a non zero minimum, 
and then increases continuously to 00. The field x decreases continuously to 
0. Correspondingly, fl decreases continuosly to flo, given in (^). 

Also, the curvature invariants all remain finite for —00 <t< 00. Hence, 
in the present model, the evolution of the toy universe with p 7^ but 7 = | 
is singularity free. 

When the initial conditions are different, the evolution can again be anal- 
ysed along similar lines. However, the main result that the evolution is 
singularity free remains unchanged. 

5 c. 7 = 1, t>U 

Initially, we have A(tj) > 0, x{U) > 0, and flfa) > 3. Then, e = 1 
necessarily. For t > U, the evolution proceeds as described in section 4a. 
Both e A and x increase. Hence, fl increases. 

The constant c > in (|15|). However, (1 — 37) < and, hence, a(t) < 0. 
Therefore, {o~(t) + c) can become negative for x < 1, or remain positive as 
X — > 1. In the case where (u(t) + c) becomes negative for x < 1, it must pass 
through a zero where x vanishes. This critical point is a maximum of x(t) 
beyond which we have A > and x < 0, and Q > 3. 

Consider the case where (cr(£) + c) remains positive as x ~ > 1- Further- 
more, as described in section 4, it can be taken to be non infinitesimal by 
choosing the model dependent constant fli in (0) large enough. The dynam- 
ics near % = 1, equivalently 4> = in our normalisation, is clear in terms of 
(f> and tp((p)- It follows from the analysis of section 4a that (p crosses and 
continues to increase. Then, ifi(<fi) and, hence, fl decrease. We then have 
A > 0, but x < 0. 

Thus, in both of the above cases, we now have A > and x < 0. The 
evolution then proceeds as described in section 4b. e A increases whereas x 
and, hence, fl decrease. As t — > 00, e A and x evolve as given in (|35D and 
©• 
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In particular, it can be seen that the quantities in (|53| ) are all finite for 
ti < t < oo, implying that all the curvature invariants are finite. Hence, the 
evolution is singularity free. 

t < U, equivalent ly (— t) > (-U) 

To describe the evolution for t < ti, we reverse the direction of time. 
Then, in terms of the reversed time variable which is also denoted as t, we 
have A(ti) < 0, x(ti) < 0, and O(tj) > 3. Then, e = — 1 necessarily. For 
t > ti, the evolution proceeds as described in section 4c. Both e A and \ 
decrease. Hence, Q decreases. Then, as shown in section 4c, there exists 
a time, say t = t m > t, t where K{t m ) = 1 in equation (|T7D, e A ( tm > > 0, 
x(t m ) > 0, and x(t m ) < 0. Note that Q(t m ) < 3. 

As shown in section 4c, K{t m ) = 1 implies that A(t m ) = 0. Hence, the 
scale factor e A reaches a minimum. For t > t m , we then have A > 0, x < 0, 
Q < 3, and e = — 1. The evolution proceeds as described in section 4d. 

It follows that a(t) = since 7 = 1. Note that the constant c < in 
(fL5"D . Thus, (er(i) + c) < c < and, hence, x < for t > t^. Therefore, x 
decreases. Eventually, x ~^ 0- Moreover, and as can be verified a posteriori, 
a{t) — > (constant) as x ^ 0. Hence, in equation (^8|), the constant c e is 
negative and non infinitesimal. Therefore, as described in section 4d, e A 
increases whereas x an d, hence, Q decrease. As t — > oo, e A and x evolve as 



given in (|4lD 



In particular, it can be seen that, for Qo ^ §, the quantities in (^) are 
all finite for U <t < oo, implying that all the curvature invariants are finite. 
Hence, the evolution is singularity free. 

Thus, the evolution of the toy universe containing only massless scalar 
field, i.e. p ^ but 7 = 1, and with the initial conditions A(ti) > 0, 
x(ti) > 0, and > 3 proceeds in the present model as follows. For t > U, 
the scale factor e A increases continuously to oo. The field x increases, reaches 
a maximum < 1, and then decreases continuously to 0. Correspondingly, Q 
increases, reaches a maximum < oo, and then decreases continuously to fl , 
given in (j?]). 

For t < ti, the scale factor e A decreases, reaches a non zero minimum, 
and then increases continuously to oo. The field x decreases continuously to 
0. Correspondingly, Q decreases continuously to Qq, given in (0). 
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Also, the curvature invariants all remain finite for — oo < t < oo. Hence, 
in the present model, the evolution of the toy universe with p ^ but 7 = 1 
is singularity free. 

When the initial conditions are different, the evolution can again be anal- 
ysed along similar lines. However, the main result that the evolution is 
singularity free remains unchanged. 

6. Evolution of observed Universe 

We now use the results of the analysis in section 4 to describe the generic 
evolution of a realistic universe, such as our observed one, in the present 
model. A realistic universe may contain different types of "matter" described 
by different values of 7, but with — 1 < 7 < 1. For example, 7 = — 1, 0, |, 
and 1 respectively for vacuum energy, dust, radiation, and massless scalar 
field. Our observed universe is certainly known to contain dust and radiation. 
However, in a variety of models, vacuum energy and/or massless scalar field 
are often assumed to be present. Our assumption that — 1 < 7 < 1 then 
covers all these cases. As follows from ([12]), "matter" with larger 7 dominates 
the evolution for smaller value of e A and vice versa. Therefore, when e A is 
decreasing we take 7 > | eventually, and when e A is increasing we take 7 < 
eventually. 

In considering the evolution of the universe, we start with an initial time 
ti, corresponding to a temperature, say > 10 16 GeV, such that GUT symmetry 
breaking, inflation, and other (matter) model dependent phenomena may 
occur for t > ti only. Our observed universe is expanding at ti and, hence, 
A(ti) > 0. For the sake of definiteness, we take x{U) > and fl(ti) > 3, 
equivalently u>(ti) > 0, as commonly assumed. 

We first describe the evolution for t > ti taking, as initial conditions 
A{ti) > 0, x(U) > 0, and Q(ti) > 3. Then, e = 1 necessarily. To describe 
the evolution for t < ti, we reverse the direction of time and take, as initial 
conditions, A{ti) < 0, x(U) < 0, and Q(ti) > 3. Then, e = — 1 necessarily. 
The required evolution is that for t > ti in terms of the reversed time variable, 
which is also denoted as t. The results of section 4 can then be applied 
directly. 

When the initial conditions are different, the evolution can again be anal- 
ysed along similar lines as follows. However, the main result that the evolu- 
tion is singularity free remains unchanged. 
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t > u 



Initially, we have A(ti) > 0, x(t») > 0, and Q(ti) > 3. Then, e = 1 
necessarily. For t > U, the evolution proceeds as described in section 4a. 
Both e A and x increase. Hence, fl increases. Eventually, as e A increases, we 
can take 7 < 0. Then, as t — > 00, e A and x evolve as given in (|3~ll). 

In particular, it can be seen that the quantities in (|53|) are all finite for 
ti < t < 00, implying that all the curvature invariants are finite. Hence, the 
evolution is singularity free. 

Thus, in the present day universe in our model, x{today) —* 1 and 
Vt{today) -> 00. Also, as follows from (0), ^^(today) = -gf(l - x) 4q_1 -> 
0. Therefore, our model satisfies the observational constraints imposed by so- 
lar system experiments, namely Vt[today) > 2000 and -^^(today) < 0.0002 
(see @ pg. 117, 124-5, and 339). 

t < U, equivalently (— t) > (—ti) 

To describe the evolution for t < ti, we reverse the direction of time. 
Then, in terms of the reversed time variable, also denoted as t, we have 
A(ti) < 0, x(ti) < 0, and fi(£j) > 3. Then, e = — 1 necessarily. For t > ti, 
the evolution proceeds as described in section 4c. Both e A and x decrease. 
Hence, Q decreases. Then, as shown in section 4c, there exists a time, say 
t = t m > U where K(t m ) = 1 in equation ([17]), e A( * m) > 0, x(0 > 0, and 
X{t m ) < 0. Note that Q(t m ) < 3. 

As shown in section 4c, K(t m ) = 1 implies that A(t m ) = 0. Hence, the 
scale factor e A reaches a minimum. For t > t m , we then have A > 0, x < 0) 
Q < 3, and e = — 1. The evolution proceeds as described in section 4d. Now, 
however, the evolution for t > t m is complicated since the universe is known 
to contain "matter" with 7 = 0. Nevertheless, all qualitative features of its 
evolution can be obtained using the results of section 4. 

For t <; t m , A > 0, x < Oj an d (&(t) + c) < 0. Hence, e A increases and 
X decreases and, eventually, "matter" with 7 < dominates the evolution. 
Then (1 — 37) > and a(t) begins to increase. Hence (cr(t) + c), which is 
negative, also begins to increase. For 7 < 0, a(t) g rows atleast as fast as t, 
as follows from flTJp. Then eventually, as described in section 4d, e A reaches 
a maximum eventually at time, say t = t\. Also, {cr(ti) + c) < and, hence, 
x(ti) < 0. Note that this critical point of e A exists independent of the details 
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of the model, as long as "matter" with 7 < is present. Further evolution 
then proceeds as described in section 4e. 

For t > tx, we then have A < 0, x < 0, and Q < 3. Also, (a(t) + c) < 
but a(t) remains increasing. Hence, depending now on the details of the 
model, (<j(t) + c) may remain negative for all t > t\, or it may reach a zero 
and become positive. 

In the first case, (cr(t) + c) remains negative for all t > t\. Hence, x < 
and x continues to decrease. The scale factor shown in section 4e, 

reaches a minimum at time, say t = tj\t > ti, with x(^m) > non vanishing. 
For t > tM, e A increases, and we have A > 0, x < 0, and Q < 3. Further 
evolution then proceeds as described in section 4d. 

Note that this means that the scale factor increases and reaches a maxi- 
mum, then decreases and reaches a minimum, then increases and so on. The 
existence of maxima of e A is model independent as long as "matter" with 
7 < is present, which is the case for the observed universe. The minima of 
e A are all non zero for the reasons described in section 4c. Their existence de- 
pends on whether (a(t) + c) remains negative or not, but is otherwise model 
independent as long as "matter" with 7 > | is present, which is the case 
for the observed universe. The evolution can thus become repetetive and 
oscillatory. 

In the second case, (<r(t) + c) reaches a zero at time, say t = t m i > t\ and 
becomes positive. It then follows that x(^m') — and that this is a minimum 
of x- F° r t > t m >, we have A < 0, x > 0, and Q < 3. Hence, e A decreases 
and x increases. Eventually, we can take 7 > |. 

Now, (1 — 37) < and a{t) begins to decrease or remains constant. Hence 
(<j(t) + c), which is positive, also begins to decrease or remains constant. 
Thus, depending on the details of po and 7, (a(t) +c) may or may not vanish 
with x < 1- H ( a (t) + c ) vanishes with x < 1, then x vanishes. This critical 
point is a maximum of x, beyond which we have A < 0, x < 0, and e = — 1. 
Further evolution then proceeds as described in section 4c if Q > 3 and as in 
section 4e if Q < 3. The evolution can thus become repetetive and oscillatory. 

If (a(t)+c) does not vanish and remain positive for x < 1 then, eventually, 
X — > 1- It follows from equation ([15]) that fix. 2 °^ e~ 6A , whereas p is given 
by (JT2|). Note that A < and e A is decreasing. Then, in this limit, the — ^~ 



term in equation (|i~0|) dominates or, if 7 = 1, is of the order of the ^ term. 



Thus we have A < 0, x > 0, x ~ > I5 equivalently — > 00, and e = — 1. 
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Also, -p- term in equation ( IQ|) dominates or is of the order of the ^ term. 
But this is precisely the time reversed version of the evolution analysed in 
section 4b, where the dynamics is clear in terms of and ip{4>)- Applying the 
results of section 4b, it then follows that 4> will cross the value after which 

and, hence, x wiU begin to decrease. We then have A < 0, x < 0, > 3. 
The evolution then proceeds as described in section 4c. The evolution can 
thus become repetetive and oscillatory. 

Thus it is clear that depending on the details of the model, the universe 
undergoes oscillations, perhaps infinitely many. As follows from the above 
description, the oscillations can stop, if at all, only in the limit x ~ * 0- The 
solutions then are given by (|4lD - (f43"l). In particular, however, the quantities 
in (|o^|) all remain finite during the oscillations, implying that all the curvature 



invariants remain finite. They also remain finite in the solutions (£|1]) - (f43|) 
if Q < | in (0). 

Thus, it follows that if Q < ~ then the quantities in (|53"D are all finite 
for ti < t < oo, implying that all the curvature invariants are finite. Hence, 
the evolution is singularity free. 

In summary, the evolution of a realistic universe, such as our observed 
one, with the initial conditions A(U) > 0, x(U) > 0, and Q(ti) > 3 proceeds 
in the present model as follows. For t > ti, the scale factor e A increases 
continuously to oo. The field x increases continuously to 1. Correspondingly, 
Q increases continuously to oo. 

For t < ti, the scale factor e A decreases and reaches a non zero mini- 
mum. It then increases and reaches a maximum, then decreases and reaches 
a minimum, then increases and so on, perhaps ad infintum. The oscillations 
can stop, if at all, only in the limit % ~~ * 0- The solutions then are given by 
( |4lD - fl4~3|) . The field x ma Y; depending on the details of model, continuously 
decrease to 0, or undergo oscillations, its maxima always being < 1. 

Also, the curvature invariants all remain finite for — oo <t< oo. Hence, 
the evolution, although complicated and model dependent in details, is com- 
pletely singularity free. Thus we have that a homogeneous isotropic universe, 
such as our observed one, evolves in the present model with no big bang or 
any other singularity. The time continues indefinitely into the past and the 
future, without encountering any singularity. 

When the initial conditions are different, the evolution can again be anal- 
ysed along similar lines. However, the main result that the evolution is 
singularity free remains unchanged. 
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7. Generalisations of the model 

The analysis and the results of sections 4-6 are generic and valid for any 
function Q satisfying only the constraints (0). However, these constraints 
can be relaxed further thus generalising the model, but still maintaining the 
evolution singularity free. We now study these generalisations. 

We first summarise the essential features of the constraints (0) on Q(x) by 
noting the following points. First, as evident from the results of the previous 
sections, any singularity is likely to arise, if at all, only in the limit x ~ * 0- 
It then follows from the relevent solutions Q4"T| ) - ( pE3| ) that this potential 
singularity is avoided by the constraint Q < ~ in (0). 

Second, we used only the sufficiency condition given in the Appendix to 
determine the absence of singularity. Hence, the requirement of finiteness of 
P.Vn>l. 

Third, there is no singularity in the limit ^ — > 1, although Q —>■ oo. The 
evolution of the present day universe corresponds to this limit. Also, the 
behaviour of fi(x) i n this limit, in particular the constraint on a given in (|7]) 
which follows naturally from (|6]), ensures that our model satisfies the obser- 
vational constraints imposed by solar system experiments. Note, however, 
that the observational constraints are satisfied for any a > 4. 

Fourth, the monotonicity of £l(x) f° r < x < 1 is n °t required either 
to avoid the singularity or to conform with the observations. It was invoked 
only to keep the analysis simple and transparent. Without this requirement, 
the evolution would be more complicated, its history possibly bearing the 
imprints of the ups and downs of Q{x)- The analysis in section 4 would then 
involve further subclasses not directly relevent for the issue of singularity. 

It is now clear how the present model, equivalently the constraints (^) on 
Q(x), can be generalised further. In the generalised model also we take x t° 
have a finite upper bound x < Xmax < oo which was naturally the case in 
the previous model. Then, as in section 2, the factor Xmax can be absorbed 
into the Newton's constant Gn and the range of x can set to be < x < 1. 
The generalised model is now given by the function fl(x) > satisfying the 
following more general constraints: 

1 

< - 

~ 3 

= finite Vn>l, 0<x<l 



fi(0) 
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— > oo at x = 1 only 
limfi = ^(l-x)^, «>~, (44) 

where f2i > is a constant. The function ft(x) is arbitrary otherwise. 

Thus, in the generalised model, fl(x) i n the limit x ~ > need not nec- 
essarily be a constant, but can be a function of X- F° r example, one may 
have fi(x) = Qo'X 13 as X ~ *• where /3 > and f2 ' > are constants. Such 
functions do arise in the variable mass theories of || also. Note that the 
position of the pole, namely x = 1, is n °t a constraint but only a choice 
of normalisation of X- However, the order of the pole a can now be, for 
example, > 1. 

We have taken ^4r, \fn > 1 for the sake of simplicity. This requirement 
can perhaps be relaxed if one uses only the appropraite necessary conditions 
to determine the absence of singularity instead of the sufficiency conditions 
given in the Appendix. However, we will not pursue it further in this paper. 

The evolution of the universe in the generalised model can be analysed 
along the same lines as in sections 4-6. It is clear that the only important 
limits are x ~ > an d, perhaps, x — ► 1. Hence, we now analysis only these 
two limits in the generalised model. To be definite, we consider the following 
example: 

lim ft = fio'X' 3 , P > 

limQ = Qi(l-x)~ 2a , a>l, (45) 

where ftc > and fti > are constants. Note that the analysis in section 
4 applies directly with no modification in the limit x — » 1 if t < en < 1. 
Hence, we have taken a > 1 in (£5J). Using @, the constraints (|45l ) can be 
expressed equivalently in terms of and ip{4>) as follows: 

2 

lim if} = -ln(0-0 o ) 
lim i/j = ipo e ~^ ; if a = 1 

= -00 , if a > 1 , (46) 

where <fto and are constants. The range of ip — lnx is given by — oo < ip < 
0. Hence, the constant i/jq < in (pj). Also, the limit x ~~ *• corresponds 



31 



to ip — > — oo and, hence, to — > 0o whereas the limit x 1 corresponds to 
■0 — > and, hence, to — > oo. 

X^O 

In this limit, the analysis is similar to that given in section 4d. Under 
the conditions of section 4b, equation ( |T7| ) can be solved in the limit x ~~ 
relating x an d ej4; 

e A = (constant) x _5 . (47) 

Hence, e A — > 00 as x ~^ 0. Substituting ( fT7|) in equation (p~5|) then yields the 
unique solution in the limit x — > 0: For 0^1, 

e A = e A ° (t - sign(m)t)"^ , X = Xo (*o - sign(m)t) m , (48) 



where t is measured in appropriate units, Aq, xch and to > U are constants 
J- 

e A = e Ao e -^ 5 x = Xoe c(t-t ) ? (49) 



and m = For /? = 1, 



where c e < is a constant. Note, as pointed out in |L5[ and as discussed in 
section 3, that when x and e A satisfies equation fl4"7|) the field equation (13) 
is not guaranteed to follow from (p~Q|) - ([T2]). Nevertheless, the above solution 
satisfies equation ( |i~3| ) also, as checked explicitly. 

Let /3 > 1. Then, m > 0. Since x — 0, it follows from (|48| ) that 
t ^ to > ti, implying that x vanishes at a finite time £o- I n this limit, the 
scale factor e A — >• 00 as follows from (48). The quantities in (|53D, for example 



^, also diverge implying that the curvature invariants, including the Ricci 
scalar, diverge. Thus, for (3 > 1, there is a singularity at a finite time to- This 
implies, for example, that the variables mass theories of which exhibit a 
minimum of e A can nevertheless have singularities because (3 can be > 1 in 
these theories. 

Let (3 < 1. When (3 < 1, m < 0. Since x — ► 0, it follows from (g8p that 
£ — > 00. In this limit, the scale factor e A — > 00 as follows now from (f48"|). The 
same result holds when (3=1 also, as follows from (f4"5p. 

If p 7^ then the consistency of the above solutions requires that 7 > |. 
This is analogous to the condition 7 > required in section 4d. 

It can now be seen, for (3 < 1, that the quantities in (^) are all fi- 
nite implying that all the curvature invariants are finite. Thus, there is no 
singularity in the above solutions when (3 < 1. 
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In this limit, the analysis is similar to that given in section 4a. For p ^ 
and 7^1, the solution is given uniquely by 



e A = e A H3(i+-y) 
X = 1 - Xoe kT , if a = 1 

= l-xo(l-fcT)" 5 ^ , if «>1, (50) 

where x > and A; > are constants, t is measured in appropriate units, 

1-7 

and we have defined T = t 1+ t which — > in the limit £ — ► oo. The above 
solution then implies that e A —* oo and % — ► (1 — % )- The precise value 
of xo is model dependent. In particular, it also depends on Qi and can be 
made as small as necessary by choosing a sufficiently large Q\. Equivalently, 
as £ — > oo, Q can be made as large as necessary. 

This will imply, as in section 6, that in the present day universe in the gen- 
eralised model, x(today) — > 1 — Xo, Q(today) — > Q^q 20 , and ^^(today) — > 
— fpXo a l - Therefore, by choosing Oi sufficiently large one can ensure that 
the generalised model also satisfies the observational constraints imposed 
by solar system experiments, namely VL(today) > 2000 and -^^(today) < 
0.0002 (see § pg. 117, 124-5, and 339). 

For p = 0, the solution is obtained easily in terms of <p and if) ((f)), and is 
given by 

e A = e A H \ 

= (constant) , (51) 

where Aq is a constant. ip(4>) and, hence, x = are then obtained from 
(f46|). Note that, as in section 4a, these solutions also apply to the case when 
p 7^ 0, but 7 = 1. As t — ► oo in the above solution, we have e A — ► oo and 



4> — ► oo. Thus, as follows from fl46|), if) — > 0. Hence, x — > 1 and — ► oo. 

It can be seen that the quantities in (|53| ) are all finite implying that all 
the curvature invariants are finite. Thus, there is no singularity in the above 
solutions. 

Now the evolution of the universe in the generalised model can be anal- 
ysed as in sections 4-6. For the three toy universes of section 5, the evolution 
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proceeds as described in section 5. For a realistic universe, such as our ob- 
served one, the evolution is complicated but proceeds as described in section 
6. In particular, it follows that the evolution of the universe is singularity 
free in the generalised model when the function fl(x) satisfies the constraints 



(45|). Although not proved here, the evolution is likely to be singularity free 
when fl(x) satisfies the constraints (f44|) only. 

Before concluding this section, we would like to point out an interesting 
aspect of the generalised models, where fi(x) ~~ > 0, equivalently uj — > — |, in 
the limit \ — ► 0. When cu — — |, the graviton-dilaton part of the action in 
(01) acquires a conformal symmetry under which 

g^^e^g^ and X - e^'x , (52) 



where rj(x) is an arbitrary function of space time coordinates |16|]. In general 



the matter action 5m breaks this symmetry. However, it can be seen from 
the solutions presented in section 4 and the present one that precisely in 
the limit X — * do the matter terms become negligible compared to the 
dilatonic terms. This suggests that the matter action Sm may be neglected 
and that the conformal symmetry of the action (HJ) may become exact as 
X — > 0. A possible implication of the appearance of this conformal symmetry 
is mentioned in the conclusion. 



8. Summary and Conclusion 

To summarise, we considered the evolution of a flat homogeneous isotropic 
universe in a graviton-dilaton model. The model is specified by the function 
fi(x) which satisfies the constraints given in (|7]), or more generally fl44|), but 
is arbitrary otherwise. We studied the evolution of universe in this model 
with the main purpose of determining whether it is singular or not. 

Assuming generic initial conditions where A(tj) and x(ti) are finite and 
non infinitesimal and < oo, we presented a general analysis which is 

valid for any matter and for any arbitrary function O(x), and which en- 
ables one to analyse the evolution and obtain its generic features even in the 
absence of explicit solutions. 

We illustrated our method by applying it to describe the evolution of three 
toy universes. We showed, in particular, that their evolution in the present 
model is singularity free. We then described the evolution of a universe, such 
as our observed one. We showed, in particular, that its evolution, although 
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rich and complicated, is singularity free. This is a generic result valid for any 
fi(x) satisfying only the constraints given in ([/]) or ([44]). 

An important question to ask now, from our perspective, is whether a 
function f2(x) as required in the present model can arise from a fundamental 
theory. We now discuss critically a few such possibilities in string theory. 

In the string theoretic context, Einstein metric formulation given in (|l|) 
is more natural where (f> is the dilaton and x(</>) the arbitrary function. The 
strong coupling limit in string theory, relevant for singularities, corresponds 
to diverging effective Newton's constant (= gp) and, thus, to the limit 

x(</>) — > 0. In this limit ln% — > -^== in our model (see equation (|7|) and 

(§)• 

(1) Note that, in string theory, ln(x{4>)) can be thought of as Kahler 
potential for (ft, which is expected to be modified at strong coupling by non 
perturbative effects But, such modifications are usually exponential in 



nature and, hence, are unlikely to lead to a x{4>) as required here. 

(2) However, at strong coupling, <fi here may not be the stringy dilaton 
directly, but instead a combination of the stringy dilaton and other compact- 
ification dependent moduli. In that case, after writing the relevant effective 
action as given in ([!]), the Kahler potential ln% for this 'effective dilaton' (j) 
may well turn out to be of the required form. 

(3) Another, perhaps more promising, possibility is the following which 
arises when Q = | (see equation (0) and (fj)). This case corresponds to a five 



dimensional space time compactified to four dimensions on a circle In 
recent developments in string theory at strong coupling, analogous phenom- 
ena relating a d-dimensional theory and a {d + l)-dimensional theory on a 
circle are found to occur: Using S-duality symmetries, Witten has discovered 
that the ten dimensional string at strong coupling is an eleven dimensional 
(M-)theory compactified on a circle [|19|, |]]. A similar phenomenon, that the 
three dimensional string at strong coupling is a four dimensional theory com- 
pactified on a circle, is at the heart of Witten's novel proposal for solving 
the cosmological constant problem p0| ; for its possible stringy realisation see 



21]. A similar phenomenon in four/five dimensions, if exists, may possibly 
lead to a x{4>) as required here. 

(4) As discussed in section 7, the graviton-dilaton action in (f|) acquires 
a conformal symmetry in the limit x ~~ * where fi(x) — > 0, equivalently 
00 (x) ~ * ~ § ■ Note that this is also the limit where singularities can poten- 
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tially arise, but are avoided in the present model. The apperance of this 
conformal symmetry is, perhaps, a hint of a connection between the present 
model and string theory, known to lead generically to conformal field theories, 
but which must now be in the strong coupling limit since x ~~ * 0- 

Admittedly, at present, these are plausibility arguments only. Neverthe- 
less, given the elegant way the present model resolves the big bang singularity 
it is worthwhile to derive it, perhaps along the above lines, from a funda- 
mental theory such as string theory. 

On the other hand, however, such a model, even if not derivable from 
string theory, may stand on its own as an interesting model for singularity free 
cosmology. One may then study its implications in other cosmological and 
astrophysical contexts. There are many issues that can be studied further. 
We mention some of them here. 

(1) One can study the evolution of a more complicated universe, e.g. an 
inhomgeneous anisotropic universe, and determine whether it is singularity 
free or not. 

(2) In section 6 we found that the evolution of a realistic universe can be 
oscillatory. Such early universe oscillations have significant implications and 
may also lead to observable predictions. It is therefore important to study 
this issue in more detail. 

(3) In the astrophysical context, one may study the evolution and/or col- 
lapse of stars. In fact, in |TTJ where the present model was originally derived 
in a different context, we found that our model suggests a novel scenario for 
stellar collapse in which a black hole is unlikely to form. Such a phenomenon, 
if established rigorously, is likely to have far reaching consequences. 

(4) The present model is characterised by one arbitrary function fl(x), 
required to satisfy the constraints given in ([?]) or flip). These constraints, 
however, do not fix fl(x) uniquely. The consequent arbitrariness in the model 
greatly diminishes its predictive power. It is therefore desirable to find the 
criteria which can fix fi(x) uniquely. 

In this regard, note that the present model naturally incorporates, as 
described in [|12||, an ingredient crucial for the success of hyperextended in- 



flation ||. Therefore, with a view to fix f2(x) uniquely, one may require the 
present model to lead to a successful inflation also. However, our preliminary 
calculations |22| indicate that this requirement is not strong enough to fix 
uniquely. Hence, more criteria are needed. 
Considering the simplicity of the present model and the novel conse- 
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quences that follow from it generically, we believe that its further study 
is fruitful and is likely to lead to interesting phenomena. 

Acknowledgement: We thank T. R. Seshadri for discussions. 

APPENDIX 

Finiteness of Curvature Invariants 

In this Appendix, we derive a sufficient condition for all curvature invari- 
ants to be finite. Curvature invariants are constructed using metric tensor, 
Riemann tensor, and covariant derivatives. When the metric is diagonal, 
every term in any curvature invariant can be grouped into three types of 

factors: (A) y/g^g vv g xx g TT R^ Xt , (B) V9^9 uu 9xx 1% and (C) ^ <9 M 
acting multiply on (A) and (B) type factors (no summation over repeated 
indices) . 

Using the metric given by (|Sp and evaluating the above factors explicitly, 
one obtains that (A) and (B) type factors are functions of A, A, and e~ A . 
The action of (C) produces extra time derivatives. Thus, it follows that the 
curvature invariants are functions of e~ A and ^4-, n > 1. 

By repeated use of equations (|il3l)-(|i~2j), C^- for any n > 1 and, hence, any 
curvature invariant, can be expressed in terms of the following quantities: 

e~ A - *, and ^ — (*Y , Vn > 1 . (53) 

X X^' X X 2 ' dx n \XJ 

The resultant expressions are finite if the above quantities are finite. The 
algebra is straightforward and, hence, we omit the details here. 

Hence, a sufficient condition for all curvature invariants to be finite is 
that the quantities in j\53j ) be all finite. 
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